Constraint Satisfaction over Generalized Staircase Constraints 
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Abstract 

One of the key research interests in the area of Constraint 
Satisfaction Problem (CSP) is to identify tractable classes 
of constraints and develop efficient solutions for them. 
In this paper, we introduce generalized staircase (GS) 
constraints which is an important generalization of one 
such tractable class found in the literature, namely, stair- 
case constraints. GS constraints are of two kinds, down 
staircase (DS) and up staircase (US). We first examine 
several properties of GS constraints, and then show that 
arc consistency is sufficient to determine a solution to 
a CSP over DS constraints. Further, we propose an 
optimal 0(cd) time and space algorithm to compute arc 
consistency for GS constraints where c is the number of 
constraints and d is the size of the largest domain. Next, 
observing that arc consistency is not necessary for solv- 
ing a DSCSP, we propose a more efficient algorithm for 
solving it. With regard to US constraints, arc consistency 
is not known to be sufficient to determine a solution, and 
therefore, methods such as path consistency or variable 
elimination are required. Since arc consistency acts as 
a subroutine for these existing methods, replacing it by 
our optimal O(cd) arc consistency algorithm produces a 
more efficient method for solving a USCSP. 

Keywords: Constraint Satisfaction Problem (CSP), Stair- 
case Constraint, Monotone Constraint, CRC Constraint, 
Generalized Staircase Constraint 



1 Introduction 

Constraint Satisfaction Problem (CSP) is one of the key 
techniques in artificial intelligence that offers a simple for- 
mal framework to represent and solve several problems in 




5 = Staircase 

MacC= Max-Closed 

MiC- Min-Closed 

CRC= Connected Row Convex 

US= Up Staircase 

DS = Down Staircase = MaC n MiC 

GS= Generalized Staircase=D5u US 



Figure 1: The space of generalized staircase constraints. 

temporal reasoning, job scheduling, pattern matching and 
natural language processing 11281 . Interestingly, in many 
such problems, the underlying constraints are of a spe- 
cific type. For this reason, one of the key research inter- 
ests in the area of CSP is to identify tractable classes of 
constraints and develop efficient solutions for them. In 
this paper, we introduce generalized staircase (GS) con- 
straints which is an important generalization of one such 
tractable class introduced in (TT), namely, staircase con- 
straints. 

Figurefllshows the relationships of GS constraints with 
other constraint classes. The class of GS constraints is 
the union of two constraint classes, namely, down stair- 
case (DS) and up staircase (US). While GS constraints 
strictly generalizes the class of staircase constraints (S, 
shown as a pie-shape), they are a strict subclass of con- 
nected row convex (CRC) constraints ATI . Although re- 
lated constraint classes such as CRC, max-closed, min- 
closed and monotone constraints have been well studied 
in the literature, to the best of our knowledge, the specific 
class of GS constraints have not been investigated earlier. 
GS constraints are interesting to study as many temporal 
constraints involving bounded intervals that arise in tem- 
poral reasoning and temporal databases lfl3l [191 [3T1 l32l 
are often GS constraints. Due to their monotonic and con- 
vex structure, they admit significantly simpler and faster 
solutions than CRC constraints. 
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Figure 2: Types of generalized staircase (GS) constraints. 

A GS constraint can be either a down staircase (DS) 
or an up staircase (US) constraint. A DS (respectively, 
US) constraint is a CRC constraint [11] such that if [p, q] 
and \p' , q'] are images of consecutive rows in its constraint 
matrix, then p < p' A q < q' (respectively, p > p' A 
q > q'). Referring to the constraint matrices shown in 
Figure [2] a DS (respectively, US) constraint has images of 
consecutive rows shifted towards right (respectively, left) 
as we move down the matrix. 

1.1 Motivation 

Consider the following simplified problem of identifying 
cyclones at a given location from weather records. Sup- 
pose a cyclonic instance is characterized by the follow- 
ing set of events: hourly rainfall exceeding a units (event 
A), hourly mean wind speed exceeding /3 units (event 
B), and clockwise wind direction (event C), such that 
the following constraints hold: (i) —3 < A — B < 1, 
(ii) -2 < B - C < 2, (iii) -2 < C - A < 3 where X 
denotes the time point of event X € {A, B, C}. These 
constraints essentially put a bound on the proximity of 
two different events irrespective of their order of occur- 
rence. For example, constraint (i) indicates that any A 
event should start in the neighborhood of B given by 
[B.s-3,B.s + l]. 

The weather database contains a list of instances of 
events A, B and C sorted on x where x denotes an in- 
stance of event X. The goal is to determine tuples of the 
form (a, b, c) where a, b, c are instances of events A, B 
and C respectively, such that all the above constraints 
are satisfied. Modeling it as a CSP, we consider vari- 
ables Xa, Xb, Xc for events A, B, C respectively, each 
of which have domains as the time points of respective 
event instances, and the constraints as shown in the con- 
straint network in Figure [3] Considering the constraint 
matrix for the first constraint, if a^ , a*' € A have images 
[6j,6fc] and [&,■',&*/] respectively in B (i.e., bj to bk sat- 
isfy the constraint for at, etc.), and i < i', then j < j' and 
k < k' (the subscripts denote the indices of the instances 
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Figure 3: Constraints for the example problem of cyclone 
identification. 

in the sorted order). This is the down staircase property as 
shown in Figure |2|a). Similarly, the other constraints (ii) 
and (iii) are also DS constraints. 

Problems such as the above are common in complex 
event processing applications ll]|4]|5]|6]|7l[8l[l9]|23]|24] 

Hi|2aim|3U|3U[33][3!|35][37) where the g0al is t0 
efficiently detect occurrences of complex events which 

are usually represented as patterns of events sharing some 
temporal relationship with each other. Examples of com- 
plex events include state changes in business and indus- 
trial processes, problems in enterprise systems and state 
changes in the environment. However, in most of the ex- 
isting literature, due to the simplicity of the patterns, they 
are represented using constructs similar to regular expres- 
sions and the detection of complex events is realized by 

using finite state automata m|5]00[U|23j|24l|26l[3l][37). 
We observe that while regular expressions can capture se- 
quential (or close to sequential) patterns quite well, they 
are inadequate to capture arbitrary pattern of temporal 
events sharing temporal relationships between pairs of 
events. Hence, we consider a CSP based approach where 
the event pattern is represented through a constraint net- 
work. In addition, we observe that the constraints in- 
volved in such constraint networks are usually DS con- 
straints, as in the above example. 

Another application of DSCSP is scheduling (or time- 
tabling) of resources (or facilities) based on their avail- 
ability. There is a variable for each resource whose do- 
main consists of the times of availability of the given re- 
source. The constraints involved are usually unary or bi- 
nary. The unary constraints are typically simple to pro- 
cess, and thus, domain values that do not satisfy those are 
easy to eliminate. The binary constraints involve bounded 
temporal intervals, which can be expressed as DS con- 
straints. Solving this constraint network would offer a 
feasible schedule for the resources. Such constraint-based 
timetabling have been considered in (9] [Kjl [TTl . 

1.2 Contributions 

In this paper, we make the following contributions: 
1 . We examine the structure and properties of GS con- 



straints, besides studying its relationships with the 
existing constraint classes. 

2. We show that, for DS constraints, arc consistency is 
sufficient to determine a solution. Further, we present 
an optimal algorithm, called ACiDS, for computing 
arc consistency for GS constraints in 0(cd) time and 
space, where there are n variables, c constraints and 
the largest domain size is d. 

3. Observing that arc consistency is not necessary for 
solving a DSCSP, we propose a more efficient algo- 
rithm, DSCSP Solver, for the same. 

4. With regard to US constraints, arc consistency is 
not sufficient to determine a solution, and therefore, 
methods such as path consistency or variable elim- 
ination are required. Since arc consistency acts as 
a subroutine for these existing methods, replacing 
any known arc consistency algorithm by our optimal 
0{cd) algorithm produces a more efficient method 
for solving a USCSP. 

1.3 Organization 

The rest of the paper is organized as follows. Section [2] 
covers the preliminary concepts of constraint satisfaction 
and formally defines the class of GS constraints. Sectionljj 
outlines the related work in this area. We next discuss 
the relationships of GS constraints with the existing con- 
straint classes and its closure properties in Section |4] In 
Section 15] we first show that arc consistency is sufficient 
to solve a DSCSP and later present an efficient arc con- 
sistency based algorithm, ACiDS, to solve a DSCSP. We 
also discuss how the ACiDS algorithm is useful to solve a 
USCSP, although arc consistency is not sufficient to solve 
it. In Section [6] we present a more efficient algorithm, 
DSCSP Solver, to solve a DSCSP. Finally, we conclude in 
Section [7] 



2 Preliminaries 

2.1 CSP over GS Constraints 

A finite binary constraint satisfaction problem (CSP) 
V = (X,T>,C) is defined as a set of n variables 
X = {Xi, X2, ■ ■ ■ ,X n }, a set of domains T> = 
{Di, £>2, ■ ■ • , D n } where Di is the finite set of possible 
values that can be assigned to Xi, and a set of c con- 
straints C = {Cij} where Cy is a binary constraint in- 
volving Xi and Xj. We assume that the domain of each 
variable consists of at most d distinct integers sorted in 



an increasing order. A constraint Cij on the ordered set 
of variables (Xi,Xj) specifies the admissible combina- 
tions of values for Xi and Xj. If a pair of values (uj, Vj) 
(vi £ Di, Vj £ Dj) satisfy the constraint Cij, then we de- 
note it as (vi, Vj) £ Cij. For the sake of simplicity, we as- 
sume that there exists at most one constraint between any 
pair of variables. Given a constraint Cij, the constraint 
Cji refers to a transposition of Cij, i.e., if (vi,Vj) £ Cij, 
then (vj, v^ £ Cji. Further we assume that all constraints 
are binary, i.e., they involve only two variables. It was 
stated in |fT6l that any r-ary CSP for any fixed r > 1 can 
be converted to an equivalent binary CSP [25 1. Therefore, 
it is reasonable to consider CSPs with only binary con- 
straints. A binary constraint can always be represented by 
a constraint matrix [21 1, showing the admissible combi- 
nations of values. 

A CSP V asks whether there exists an assignment of 
the variables X such that all the constraints are satisfied. 
If an assignment X t — Vi , Vi is a solution, then Vi is said 
to be a member of a solution. 

Consider the set of vertices Vp = {Vi, V2, ■ ■ ■ ,V n } 
where Vi corresponds to variable Xi. A pair of variables 
(Xi , Xj ) that has a constraint Cij is referred to by the arc 
(i.j) which represents an edge between Vi and Vj. The 
set of all arcs is denoted by arc(V). The graph induced 
by the vertex set Vp and the edge set arc(V) is referred 
to as the constraint network of V and is denoted by Np, 

Following [11|, a constraint CV, is row-convex if and 
only if in each row of the matrix representation of Cij, 
all the ones are consecutive, i.e., no two ones within 
a single row are separated by a zero in that same 
row. The reduced form of a constraint Cij, denoted by 
C*j, is obtained by removing all the empty rows and 
columns in its matrix representation. The domain of Xj 
through the constraint CV,, denoted by Di(Cij), is the 
set {vi £ Di\3vj : («j, Vj) £ Cij}. Suppose Cij is a row- 
convex constraint and v% £ Di(Cij). The image of Vi in 
dj is the set {vj £ Dj | (uj, Vj) £ Cij}. Due to row con- 
vexity of C^, this set can be represented as an interval 
[wi, w m ] (over the domain Dj(Cji)) and we denote w\ 
and w m by min(Cij,Vj) and max(Cij, Vj) respectively. 
We denote by succ(vj, Dj(Cji)) and pred(vj,Dj(Cji)) 
the successor and the predecessor of Vj in Dj(Cji) 
respectively. For ease of notation, we simply use 
min(vj), max(vj), succ(vj) and pred(vj) when there is 
no ambiguity on the underlying domain. 

Referring to the example in Figure [2] a), min(a$) = 
b3,max(a3) = &4,succ(&3) = 64 and pred(&3) = hi- 
A row-convex constraint CV; is connected row-convex 
(CRC) if and only if the images [a, b] and [a' , b'] of two 
consecutive rows in C*- are such that a' < succ(b) A b' > 
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Figure 4: (a, /3)-monotone constraints. 

pred(a). Following [30|, a constraint CV, is said to be 
monotone if and only if (v,w) G C,-j implies {v'^v'A £ 
C,j for all values u,,^ € -Di, Uj,u'- £ -Dj such that 
w- > Uj and i>' < Vj, i.e., referring to the constraint ma- 
trix of Cij, if («j, u,-) € Cy, then all cells to the left of 
(vi , Vj ) on the same row are in Cij and all cells below 
(vi, Vj) on the same column are in Cij. However, a com- 
plementary definition of monotone constraints was pro- 
posed in |fl4l. where a constraint Cij is said to be mono- 
tone if and only if (vi, Vj) £ Cy implies (v[, v'A £ CV, 
for all values Vi,v[ € Di, Vj,v'j £ Dj such that v\ < v t 
and v'j > Vj, i.e., referring to the constraint matrix of Cij, 
if (vi,Vj) £ Cij, then all cells to the right of (vi, Vj) on 
the same row are in Cij and all cells above (vi, Vj) on the 
same column are in CV, . An example of each of them is 
shown in Figure Wa) and Figure Wb). Now we consider 
a generalization of monotone constraints, introduced in 
ifTTl . as following. 

Definition 1. Let ^ and y be total orderings on Di and 
Dj, respectively. A constraint Cij is (<,>)-monotone 
if and only if (vi,Vj) £ Cij implies (v'^v'j) £ Cij for 
all values v i ,v' i £ Di, Vj,v'a £ Dj such that v[ < v^ 
and v'j >z Vj. Staircase constraints are (a, /3)-monotone 
constraints where a, f3 £ {<, >}. 

Figure HI shows examples of staircase constraints. 

In this paper, we propose generalized staircase (GS) 
constraints as a generalization of staircase constraints. GS 
constraints are of two types, down staircase (DS) and up 
staircase (US), defined as below. 

Definition 2. A constraint Cij is a down staircase (DS) 
constraint if and only if it is row convex and for any 



u, v £ Di such that v = succ (u), the following condi- 
tions (DS property) hold: min(Cij,u) < min(Cij,v) 
and max (Cij, u) < max [Cij, v). 

Definition 3. A constraint Cij is an up staircase (US) 
constraint if and only if it is row convex and for any 
u,v £ Di such that v = succ(u), the following condi- 
tions (US property) hold: min(Cij,u) > min(Cij,v) 
and max (Cij, u) > max (Cij, v). 

Figure [2] shows examples of both the constraints. 

Following the above definitions, it is clear that this 
class of GS constraints is a strict subclass of CRC con- 
straints. Further, (<, >) -monotone and (>, <) -monotone 
constraints are strict subclasses of DS constraints while 
(<j<) -monotone and (>,>) -mono tone are strict sub- 
classes of US constraints (see Figure [TJ. A binary CSP 
with all its constraints as DS (respectively US) constraints 
is referred to as a DSCSP (respectively USCSP). 

2.2 Arc Consistency and Path Consistency 

Following |Q3), an arc(i,j) is arc consistent if and only 
if Vvi £ Di , 3vj £ Dj such that (vi , Vj ) £ Cij (vj is said 
to support Vi). For arc consistency, we assume that con- 
straint Cij £ C if and only if Cji £ C, i.e., arc(i,j) and 
arc(j, i) are treated distinct, and that both are in arc(V). 
If each arc(i,j) £ arc(V) is arc consistent, then the con- 
straint network Af-p is said to be arc consistent. 

Following [11], a constraint network Af-p is path con- 
sistent if and only if for every triple (Xi , Xk , Xj ) of vari- 
ables such that arc(i, k), arc(k,j),arc(i,j) £ arc(V), 
for every Vi £ Di and Vj £ Dj such that (v i} Vj) £ Cy, 
there exists Vk £ Dk such that (vi,Vk) £ Cik and 
(v kl Vj) £ C k j- 

Two constraint networks are equivalent if they have the 
same set of solutions. Given a constraint network, the 
task of an arc consistency (respectively path consistency) 
algorithm is to generate an equivalent constraint network 
that is arc consistent (respectively path consistent). 

3 Related Work 

Although CSPs are NP-hard, several tractable classes 
of CSPs have been identified (22). Among the no- 
table tractable classes of constraints are functional, anti- 
functional, monotone, min-closed, max-closed and con- 
nected row convex (CRC) constraints. In this paper, we 
are considering another important tractable class of gen- 
eralized staircase (GS) constraints. 

In this work, we propose an algorithm named ACiDS, 
for computing arc consistency over DS constraints. In 



this context, we note that there are several existing arc 
consistency algorithms for general arbitrary constraints. 
The most prominent ones include AC3 [18|, AC4 |20|, 
AC5 JH, AC6 and AC7 |12|. Although AC3 has 
non-optimal worst case time complexity (0(cd 3 )), it is 
still considered as the basic constraint propagation algo- 
rithm owing to its simplicity. Its successors AC4, AC5, 
AC6 and AC7 have optimal worst case time complex- 
ity of 0(cd 2 ). AC5 is a generic arc consistency algo- 
rithm that offers a framework to compute arc consistency 
for any constraint network and can be optimized for cer- 
tain classes of constraints. A refinement of AC3 in the 
form of AC2001 was proposed in |2| which has an op- 
timal worst case time complexity of 0(cd 2 ). None of 
these algorithms assume knowledge of the semantics of 
the underlying constraints. Although they agree on the 
worst case time complexity (besides AC3), their perfor- 
mance varies for different constraint classes. In addition, 
the tightness and slackness of the constraints also affect 
the performance for any given constraint class. Since arc 
consistency is sufficient to solve a DSCSP, as shown in 
Theorem|2] any DSCSP can be solved by any of these arc 
consistency algorithms (AC4, AC6, AC7 or AC2001) with 
a worst case time complexity of 0(cd 2 ). The arc consis- 
tency algorithm presented in this paper, namely ACiDS, 
has an optimal worst case time complexity of 0(cd) for 
the class of DS constraints, thus improving the 0(cd 2 ) 
bound for solving DS constraints using the existing arc 
consistency algorithms. This algorithm is a specialization 
of the AC5 algorithm. Similar specializations for several 
classes of arithmetic constraints such as functional, anti- 
functional and monotone constraints were stated in fl4*l . 
For each of these classes, AC5 determines a solution in 
O(cd) optimal time. 

GS constraints are a subclass of CRC constraints. 
While path consistency is sufficient to solve a CSP over 
CRC constraints, we show that arc consistency is suf- 
ficient to solve a CSP over DS constraints (please re- 
fer to Theorem |2j. However, it is not known that 
whether US constraints are solvable using arc consis- 
tency. For the class of CRC constraints, a path consis- 
tency based algorithm was proposed in 1 1 1 1 that runs in 
0(n?d 2 ) time. Based on variable-elimination, [36] stated 
an O (nda 2 + cd 2 ) time algorithm for CRC constraints 
where a < n is the degree of elimination of the trian- 
gulated graph of the underlying constraint network. Both 
these algorithms for CRC constraints use arc consistency 
as a subroutine which takes 0(cd 2 ) time. US constraints 
being a subclass of CRC constraints can be solved by ei- 
ther of these techniques mentioned above. Our arc con- 
sistency algorithm for GS constraints that runs in 0(cd) 



time, acting as a subroutine, thus improves the time com- 
plexity of the algorithm by 1 36 1 to O (ndo 2 + cd) (which 
is linear in d), when applied to US constraints. 

Referring to the classes of max-closed and min-closed 
constraints, [15] showed that their solutions can be com- 
puted in O (c 2 d 2 A 2 ) time where A < d is the maximum 
number of supports that a value in a domain can have. In 
this paper, we show that a constraint is a DS if and only 
if it is both max-closed and min-closed. Our DS algo- 
rithms have significantly lower time complexities for this 
restricted subclass. 



4 Properties of GS Constraints 

4.1 Relationships with other Constraint 
Classes 

In this section, we discuss the relationship of GS con- 
straints with the existing constraint classes. Figure [T] il- 
lustrates the space of GS constraints. Since some of the 
relationships of GS constraints shown in this figure have 
already been covered so far, we next explain those that 
have not yet been discussed. 

Following 03), a constraint C is max-closed 
if and only if \/(u,v),(u' ,v') £ C, implies 
(max{u,u'},max{j),t)'}) £ C. Similarly, a constraint 
C is min-closed if and only if \/{u,v),(u',v') £ C, 
implies (mm{u,u'},mm{v,v'}) £ C. The authors in 
lfT5ll claimed that if a constraint is both max-closed and 
min-closed, then it is row convex. If a constraint is either 
a (<, >) -monotone or (>, <) -monotone, then it is both 
max-closed and min-closed. We integrate these claims in 
the following result. 

Theorem 1. A constraint is down staircase if and only if 
it is both max-closed and min-closed. 

Proof. Suppose CV, is a constraint that is both max-closed 
and min-closed, and therefore, row convex |fT5l . For any 
u, v! € Di such that u' = succ(u), the images of u, v! in 
Dj are intervals, say [v, w) and [v' , w']. Since u < u', we 
claim that v < v' and w < w'. If v' < V, then using the 
min-closed property (u, v') € Cy. This contradicts the 
assumption that the image of u in Dj is [v, w}. If w' < w, 
then using the max-closed property (u',w) £ CV,, thus 
contradicting the assumption that u' has the image [v' , w'\ 
in Dj . Hence, Cy must be a down staircase constraint. 

Now, suppose Cij is a down staircase constraint. From 
the down staircase property, for any u, v! £ Di such that 
u < u', if the images of u, u' in Dj are [v, w] and [v 1 , w'), 
then v < v' and w < w'. If w < v', both max-closed and 
min-closed properties are trivially satisfied. 
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Figure 5: Examples of max-closed and min-closed con- 
straints. 

Next, consider the case v < v' < w < w', and any 
y, y' E Dj such that v < y < w and v' < y' < w'. 
If z/ < y', then the min-closed and max-closed properties 
are trivially satisfied for any pair (u,y), (u',y r ) E Cy. 
Hence, assume y > y'. However, such y' exists only for 
succ(v') < y < w. Observe that (u,y r ) E Cij for if < 
y' < y. Also (u',y) E Cy for any v' < y' < y < 
w. Hence, both min-closed and max-closed properties are 
satisfied for any pair (u, y), (v! , y') E Cij. 



□ 



Figure [T] shows the intersection space of max-closed 
and min-closed constraints as DS constraints. Now we 
consider some of the other relationships shown in the fig- 
ure. We note that there are CRC constraints that are not 
GS, but min-closed (Figure |3Ja)), or max-closed (Fig- 
ure Bib)). From Figure l5tc) and Figure BJd), we also see 
the existence of constraints that are either min-closed or 
max-closed, but not CRC. Figure |5Je) and Figure |3Jf) il- 
lustrate existence of GS constraints that are both DS and 
US. While Figure BJe) is not a staircase constraint, Fig- 
ure BJf) is staircase. Figure |4] shows that any staircase 
constraint is either max-closed, or min-closed, or both. 
This implies that any US constraint that is also a stair- 
case constraint, is either max-closed or min-closed, but 
not both. We also claim that if a US constraint is either a 
max-closed or a min-closed constraint, it must be a stair- 
case constraint. However, we note that there are US con- 
straints that are neither max-closed nor min-closed. Fig- 
ure |2jb) illustrates this fact. 

4.2 Closure Properties of GS Constraints 

Following ifTTI . we know that CRC constraints are closed 
under the following operations: transposition, intersec- 



tion and composition. Transposition of the constraint ma- 
trix implies that if (vi,Vj) E Cij, then (vj,Vi) E Cji. 
Intersection of two binary constraints Cy and Cy im- 
plies that (vi,Vj) E C^ n C'^ if and only if (vi, Vj) E 



C^ A (vi,Vj) E 



G 'iy 



Composition refers to multipli- 



cation of two constraint matrices, i.e., given two binary 
constraints Cij and Cjk, the constraint Cik — Cij x Cjk 
is such that (vt, Vk) E Cik if and only if 3vj\(vi, Vj) E 
Cij A (vj,v k ) E Cj k . 

These closure properties are necessary for path consis- 
tency which ensures tractability of CRC constraints. In 
this section, we discuss the closure properties for the sub- 
classes of GS constraints, namely DS and US constraints. 
We show that DS constraints are closed under transpo- 
sition, intersection and composition. US constraints are, 
however, closed under transposition and intersection, but 
not under composition. In fact composition of two US 
constraints result in a DS constraint. 

Lemma 1. DS constraints are closed under transposition. 

Proof. Suppose Cy is a DS constraint, and Cji is its 
transpose. 

From Theorem [T] we know that Cy is both min- 
closed and max-closed. Hence, for any two tu- 
ples (a,b),(c,d) E Cij, (min{a, c},min{6, d}) E 
Cij and (max{a,c},max{6, d}) E Cij. From the 
definition of transposition, for every (u,v) E Cij, 
(v,u) E Cji. Therefore, for any two tuples 
(b, a),(d, c) E Cji, (min{6, d}, min{a,c}) E Cji and 
(max{&, d}, max{a, c}) E Cji. Hence, Cji is both min- 
closed and max-closed. From Theorem [T] we conclude 



that Ca is a DS constraint. 



□ 



Lemma 2. DS constraints are closed under intersection. 

Proof. Suppose Cij,C^ are two DS constraints on the 
same domains Di and Dj. Assume that Aij — Cij f] C[.; 
denote the intersection of Cy and C[j. Suppose Ay is not 
a DS constraint. 

Consider any two non-empty rows of Ay, u,v E 
Di such that u < v. Suppose u, v have the re- 
spective images [a, b],[c, d] in Cy. Similarly, let u,v 
have images [a 1 , b'], [c',d'] respectively in C,'-. Note 
that min(Aij,u) — max{a, a'} and max(Aij 1 u) = 
min{6, b'}. Similarly, min(Aij,v) = max{c, c'} and 
max(Aij,v) = minjd, d'}. 

From the closure properties of CRC constraints in [ 1 1 1, 
A^ must be a CRC constraint. For Cji not to be a DS 
constraint, either of the following two cases must hold: 
(1) min(Aij,u) > mm(Aj,-,v). From the DS property 
of Cij, a < c. Similarly, a' < c! . This, however, con- 
tradicts min(Aij,u) — max{a,a'} > max{c, c'} = 



min(Aij, v). 

(2) max(Aij,u) > max(Aij,v). From the DS prop- 
erty of Cij, b < d and b' < d! . This, however, con- 
tradicts max(A i j, u) = min{6, b'} > min{d,d'} — 

max(Aij,v). D 
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Lemma 3. US constraints are closed under transposition. 

Proof. Suppose Cy is a US constraint, and Cji is its 
transpose. Suppose Cji is not a US constraint. 

Consider any two rows of matrix representation of Cji, 
say u, v E Dj such that u < v. Suppose the im- 
ages of u, v in Cji are [a, b] and [c, d] respectively where 
a, b,c,d € Dj. From the closure properties of CRC con- 
straints in ifTTI . Cji must be a CRC constraint. For Cji 
not to be a US constraint, either of the following two cases 
must hold: 

(1) a < c. From the definition of max values, v < 
max(c). Since a < c = min(v) and Cij is row con- 
vex, max(a) < v. From the above relations, it follows 
that max(a) < max(c). This, however, contradicts the 
US property of Cij . 

(2) b < d. From the definition of min values, min(b) < 
u. Since b = max(u) < d and Cij is row convex, 
u < min(d). From the above relations, it follows that 
min(b) < min(d). This, however, contradicts the US 
property of Cij . D 

Lemma 4. US constraints are closed under intersection. 

Proof. Suppose Cij , CL are two US constraints on the 
same domains D. L and Dj . Assume that Aij — Cij n C[,: 
denote the intersection of Cij and Cy . Suppose Ay is not 
a US constraint. 

Consider any two non-empty rows of Aij, u,v G D L 
such that u < v. Suppose u, v have images [a, b], [c, d] 
respectively in CV, and [a 1 , b'], [c', d'] respectively in C-.-. 
Note that min(Aij , u) — max{a, a'} and max(Aij , u) — 
min{b, b'}. Similarly, min(Aij,v) — max{c, c'} and 
max(Aij, v) = min{d, d'}. 

From the closure properties of CRC constraints in ifTTIl . 
A^ must be a CRC constraint. For Cji not to be a US 
constraint, either of the following two cases must hold: 

(1) min(Aij, u) < min(Aij,v). From the US property 
of Cij, a > c. Similarly, a' > c'. This, however, con- 
tradicts min(Aij,u) — max{o, a'} < max{c, c'} = 
min(Aij, v). 

(2) max(Aij, u) < max(Aij,v). From the US prop- 
erty of C^, b > d and b' > d! . This, however, con- 
tradicts max(Aij,u) = min{fe, b'} < min{d, d'} = 
max(Aij,v). □ 



Figure 6: Composition of two DS constraints. 

Next we show that the composition of two DS con- 
straints result in another DS constraint, as illustrated by 
an example in Figure [6] 

Lemma 5. DS constraints are closed under composition. 

Proof. Suppose Cik — A^ x Bjk is the composition of 
two DS constraints Ay and Bjk with all non-empty rows 
and columns. 

Suppose any u £ Di has image [a, b] in Aij . We claim 
that if c = min(Bjk, a) and d — max(Bji C) b), then the 
image of u in Cik is [c, d}. 

Consider any v £ Di such that v > u. By a similar 
argument, if v has an image [a' , b'] in A^, then the image 
of v in Cik is [c', d'] where c! — min(Bjk, a') and d' = 
max(Bjk, b'). 

Since v > u and Ay is a DS constraint, a' > a and 
b' > b. Using these relations and the fact that Bjk is a DS 
constraint, we infer that c < c' and d < d! . Thus, for any 
u,v G Di such that u < v, min{Cik,u) < min(Cik,v) 
and max(Cik,u) < max(Cik,v). Hence, Cik is a DS 
constraint. □ 

Next, we show that composition of two US constraints 
result in a DS constraint, as illustrated by an example in 
Figure [7] 

Lemma 6. Composition of US constraints results in a DS 
constraint. 

Proof. Suppose Co- = A^ x Bjk is the composition of 
two US constraints Ay and Bjk with all non-empty rows 
and columns. 

Suppose any u € Di has image [a, b] in Aij . We claim 
that if c = min(Bjk,b) and d = max(Bjk,a), then the 
image of u in Cik is [c, d}. 

Consider any v e D t such that v > u. By a similar 
argument, if v has an image [a 1 , b'] in A^, then the image 
of v in Cife is [c', <i'] where c' = min(Bjk, b') and <i' = 
max(Bjk,a'). 

Since v > u, and Ay is an US constraint, a! < a and 
fe' < fo. Using these relations and the fact that Bjk is 
an US constraint, we observe that c < c' and d < d' . 
Thus, for any u, v € Di such that u < v, min(Cik,u) < 
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Figure 7: Composition of two US constraints. 



there exists a smallest available support say Vj € Dj such 
that (fi,Vj) £ Cij. Similarly, since fj £ ZX-, there ex- 
ists a smallest available support say Vi £ D[ such that 
i v i, fj) G C^. From the definition of fi,fj, it follows 
that since (fi,fj) $. Cij, hence fi < Vi and /j < Vj. 
Since Cy is min-closed, (mm{fi,Vi},mm{fj,Vj}) £ 
Cij, i.e., (fi, fj) £ Cij, which is a contradiction. 

Therefore, Sf is a solution to T 5 . Similarly, we can 
show that Si is also a solution to V. □ 



min(Cik,v) and max(Cik, u) < max(Cik,v). Hence, 
Cife is a DS constraint. □ 

Owing to monotonicity and row convexity, GS con- 
straints admit much simpler and efficient algorithms to 
compute the above operations, as opposed to those for 
CRC constraints. From the arguments given in the above 
proofs, it is easy to see that each of these operations can 
be performed in 0(d) space and time for both DS and US 
constraints. Since DS constraints are solvable using arc 
consistency (as described in Section [5]), we do not need 
these closure properties. However, for US constraints, if 
we use either of the existing approaches for CRC con- 
straints (mentioned in Sectionp), we require these proper- 
ties, and hence, lies the advantage of efficient algorithms 
to compute these operations. 

5 Arc Consistency for GS Con- 
straints 



5.2 ACiDS Algorithm 

5.2.1 Algorithm Overview 

We now present the algorithm ACiDS (pseudocode in Al- 
gorithm [T] with the ArcConsf) and LocalArcConsf ) rou- 
tines in Algorithm [2]i that computes arc consistency for a 
DSCSP. This algorithm is a specialization of the generic 
AC5 algorithm lfT4l . At a high level, the algorithm first 
detects the values that are without any support: {vj £ 
Dj\$vi £ Di : (vj,Vi) £ Cji}. Then for all k such 
that Ckj £ C, the tuple (k, j, Vj) is added into a queue Q 
which essentially holds values due for constraint propaga- 
tion while pruning. After dequeuing the tuple (k,j,Vj) at 
a later stage, it checks whether Vj was the last available 
support for any Vk £ D^, and puts all such values into Q. 
The process terminates when Q becomes empty. 



5.2.2 The MIN Data Structure 



In this section, we first show that arc consistency is suffi- 
cient to determine a solution to a DSCSP, following which 
we present an algorithm, ACiDS (Arc Consistency for 
Down Staircase constraints), to compute arc consistency 
for a DSCSP. We also discuss how the same algorithm can 
be used for computing arc consistency for a USCSP. 

5.1 Sufficiency of Arc Consistency for Solv- 
ing DSCSP 

Theorem 2. Arc consistency is sufficient to determine a 
solution to a DSCSP. 

Proof. Given a DSCSP V, suppose \/i, D^D^ respec- 
tively denote the domain of Xj before and after achiev- 
ing arc consistency, and Vi, fi,k denote the first and last 
values of D' i respectively (fi < U). We claim that the 
sets Sf = {/i, / 2 , • ■ • , /«} and Si = {h,h, ■ ■ ■ Jn} are 
solutions to V . 

Suppose Sf is not a solution. Then there must exist a 
constraint Cij such that (fi,fj) $. Cij. since /,; <G D' L , 



For efficiently computing the set of values in Di for which 
Vj is the only support, we introduce a new data struc- 
ture MIN. We define MIN(vj,i) as the set of values 
in Di for which Vj is the smallest available support in Dj 
over the constraint Cij. Initially, we set MIN(vj,i) = 



{Vi\v 
V 



j — min(Cij,Vi)}. If $Vi £ Di such that 
= min(Cij,Vi), then MIN(vj,i) — <\>. In Fig- 
urega), MIN(bi,i) = [oi,o 2 ], MIN(b 3 ,i) = [0,3,0,4], 
MIN(b 5 ,i) = [05], mdMIN(b 2 ,i) = MIN(b 4 ,i) = 
4>. Due to the DS property of Cji, the set MIN(vj,i) 
can be represented as an interval. To illustrate this, 
suppose MIN(vj,i) — {vi,V3,v 6 } where Vi £ Di 
for 1 < i < 6. Since Cij is a DS constraint and 
min(v\) = min(v3) = min(ve), it implies that ^2,^4 
and V5 does not have any support in Dj over the con- 
straint Cij. Hence, we can also denote MIN in an in- 



terval form as MIN(vj 



[v\ , Ve] which essentially 



marks the boundaries of MIN(vj , i). 

We observe that since C^ is a DS constraint, for any 
Vj,v'j £ Dj such that Vj < vL if MIN(Vj,i) = [a, b] 
and MIN&Li) = [a',b'}, then b < a'. The algorithm 



Algorithm 1 ACiDS 



Algorithm 2 ArcCons and LocalArcCons 
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procedure mainCT 3 = (X,V,C)) 
Initialize^) 
for each Cy £ C do 

ArcCons(i, j, A) 
Enqueue(i, A, Q) 
Remove(A,Di) 
while Q ^ do 

Dequeue(i, j, Vj , Q) 
LocalArcCons(i , j, Vj , A) 
Enqueue(i, A, Q) 
Remove (A, D;) 

procedure Initialize 

Vu, e £>j, VCy e C, Af/JV (i)j, i)<- <£ 
for each Cy e C do 

scarii , scarij point to the first value of Di, Dj respec- 
tively 

while scarii ^A scarij ^ <f> do 
while min(scani) > scarij do 

scarij <— min(scarii) 
MIN (scarij, «)•(— Af/iV (scarij, i) U {scan;} 
scam i— succ (scarii) 

procedure Enqueue^, A, Q) 
for each u; 6 A do 

O^QU{(j, j, «<) |V (j, i) e arc (V)} 

procedure Dequeue^, j, vj , Q) 
KQj= </>then 

Q <- Q\ (i,3,vj) 

procedure Remove^, Di) 
for each dj £ C do 
if v e Di(dj) then 

pred(succ(v,Di(Cij)) ,Di(Cij)) 

■^pred(v,Di(Cij)) 
succ(pred(v,Di(Cij)) , Z)j(Cy)) 
*r-SUCc(v,Di(Cij)) 

£>i(Cy)<(-.D i (Cy)\M 



D,^D % \ {v} 



exploits this disjoint interval property of MIN sets. Dur- 
ing the course of the algorithm, the MIN sets and the 
pred, succ values get updated, while rain and max val- 
ues remain unchanged. 

5.2.3 Re-visiting ArcCons() and LocalArcCons() 

We re-design the generic subroutines ArcCons() and Lo- 
calArcConsf), stated in AC5, specifically for DSCSP in 
Algorithm [2] 

If Vi $. Di{Cij), i.e., Vi has no support in Dj 
over the constraint C\j, then we define min (Cy, fj) = 
max (Cij, Vi) = (j>. We use pointer scarii to scan the val- 
ues of Di. When scarii points to Vi, we say scarii = v. b . 
Further, scarii — <fi if and only if scarii reaches the end 
of domain Di. Similarly, we assume that succ(vi) = <f> if 
and only if Vi is the largest value in Di. 



procedure ARCCONSfi, j, A) 

A<-{ui e Di\$vj e Dj : (vi,vj) e Cy} 

procedure LocalArcCons(j, j, ty, A) 
if MIN(vj,i) = 0then 
return 

A^ <j> 

MIN (vj ,i) = [vi 1 , Vi 2 ] for some Uj, , D; 2 £ D; 
succ(?;j) 

while In'. = <f> V max (Cij , scarii) < 1/ ) h(scarii < i>; 2 ) 



scant 



do 



A-f-A U {scan;} 
scarii 4— succ(scarii) 
if scarii j^ <fi A scarii < Uj 2 then 



The algorithm begins with initialization of data struc- 
tures: queue Q to an empty queue and MIN as per the 
definition stated above. Following this, ArcCons(i,j, A) 
is called for each constraint Cij £ C. ArcCons(i,j, A) 
identifies those values V{ £ Dj such that Wj ^ Di(Cij) 
and adds them to the set A. Further, Enqueue(i, A, Q) 
enqueues tuples of the form (j,i,Vi) into Q for all 
arc(j, i) £ arc(P) and all Uj £ A. Then, the call 
to Remove(A,Di) removes all Vi £ A from D^ and 
Di{Cij) for each constraint CV,. In addition, it resets the 
pred and succ values for the preceding and succeeding 
values of u, in Di . 

The algorithm continues until Q becomes empty. 
Within the while loop (lines 7J[T0 >, a tuple (i,j, Vj) is de- 
queued from Q using Dequeue[i,j,Vj,Q) and then the 
subroutine Local ArcCons(i,j, Vj, A) is invoked. 

The heart of the ACiDS algorithm lies in the subrou- 
tine LocalArcCons (i, j, Vj, A). Firstly, it determines the 
set A = {vi £ MIN(yj,i)\succ(vj) — 4> V max(vi) < 
succ(vj)} as the set of values Vi £ Di for which Vj £ Dj 
was the only available support. Since the constraint is a 
DS, these values of A can be identified by sequentially ac- 
cessing the values in the MIN set using the succ values, 
without visiting any v ^ A. Hence, this step requires at 
most 0(A) operations. Secondly, if w' = succ(vj) ^ ^ 
it adds the set {MIN(vj,i) - A} to the set MIN(v'j,i). 
This is to ensure that those values Vi £ Di for which Vj 
was the smallest available support in Dj but not the only 
remaining support, will henceforth have v'a as the smallest 
available support. Since MIN sets are stored as disjoint 
intervals, this step can be completed in O(l) time. It is 
interesting and important to note here that if Vj is not the 
smallest support for any Vi, i.e., MIN(vj , i) = <j>, then no 
data structure needs to be updated and Local ArcConsQ 
returns immediately. 



5.2.4 Correctness of ACiDS Algorithm 

In order to establish the correctness of the ACiDS algo- 
rithm which is a specialization of AC5, we essentially 
need to show the correctness of Local ArcConsQ. In or- 
der to do so, we first show that MIN(vj , i) always retains 
the set of values in Di for which Vj is the smallest avail- 
able support in Dj over the constraint Cy . 



Lemma 7. At the end of any iteration of the set: A 
LocalArcCons(i, j,Vj, A) over a given ordered pair 
(i,j),foranyvi £ -Dj\A and Vj £ Dj, Vi £ MIN(vj,i) 



able support for Vi over the constraint Cij and now Vj is 
already pruned from Dj, Vj> is the smallest available sup- 
port for Vi at the end of the r = p th induction step. 

Thus, the induction hypothesis follows. □ 

The next lemma proves the correctness of the subrou- 
tine LocalArcConsf). 

Lemma 8. LocalArcCons(i, j, Vj, A) correctly computes 



{vi £ Di\ (Vi,Vj) £ 



Cij,Wj 



if and only if Vj is the smallest available support for V{ 
over the constraint Cij. 

Proof. We prove this lemma by applying induction on it- 
eration r. 

Consider any Vi £ D L . Suppose min(Cij,Vi) = Vj. 
Note that during initialization, w.j £ MIN(vj,i). Now, 
consider the call to LocalArcCons(i, j, Vj r , A) for the first 
time (r = 1) over the pair (i,j). If Vj 1 ^ Vj, Vi continues 
to be in MIN (vj, i) and Vj continues to be the smallest 
available support for Vi. If, however, Vj 1 — Vj, then the 
following two cases arise: 

(1) If succ(vj) — <f> or max(Cij,Vi) < succ(vj) such 
that succ(vj ) 7^ <j), then using the DS property, it implies 
that Vi lost its only support Vj and is, thus, inserted to 
A. In this case, the hypothesis is no longer applicable as 
Vi^Di\A. 

(2) If max{Cij, vi) > succ(vj) such that succ(vj) =/= <f>, 
then Vi is added to MIN{vji,i) (line 12 1 where Vj/ — 
succ(vj). By definition of succQ, Vji is the next greater 
value after Vj in Dj(Cji). Since Vj — min(Cij, vi) is al- 
ready pruned from Dj, v'j is the smallest available support 
for Vi over the constraint C^ . 

Hence, the hypothesis holds for r = 1. 

Suppose the hypothesis holds for r = p — 1, Now 
consider the call to LocalArcCons(i,j,Vj r , A) for the 
r = p th iteration over the pair (i,j). If Vj 7^ Vj, then 
following the induction hypothesis, Vi continues to be in 
MIN (vj, i) and Vj continues to be the smallest available 
support for Vi. If, however, Vj — vj, then the following 
two cases arise: 

(1) If succ{vj) — (j) or max(C'ij,Vi) < succ(vj) such 
that succ(vj ) 7^ (j>, then using the DS property, it implies 
that Vi lost its only support Vj and is, thus, inserted to 
A. In this case, the hypothesis is no longer applicable as 
Vi i A \ A. 

> succ(vj ) such that succ(vj ) 7^ <f>, 



(2) If max(Cij 

then Vi is added to MIN(vj',i) (line 12 1 where Vj 



succ(vj). By definition of succQ, Vj> is the next greater 
value after Vj in Dj(Cji). Following the induction hy- 



pothesis, for 



p — 1, since Vj was the smallest avail- 



Proof From Lemma|7| MIN(Vj,i) = <j), if and only if 
Vj is not the smallest available support for any V{ £ D^, 
Hence, pruning Vj would not make any Vi inconsistent. 
Therefore, the algorithm returns immediately. 

Now consider the case MIN(vj , i) ^ <j). Note that any 
Vi that is supported by Vj and has Vj as the only available 
support in Dj must lie in MIN(vj, i). From the pseu- 
docode (lines 5p2 1, we note the following two cases: 

(1) If succ(vj) = <p, then using Lemma[7] there exists no 
v'a larger than Vj in Dj that can potentially support any Vi. 
Hence each Vi £ MIN(vj,i) is added to A. 

(2) If succ(vj) 7^ <j>, then suppose MIN(vj,i) = 
[vi x ,..., Vi 2 ,v i3 ,...,v i4 ], such that max(v i2 ) < 
succ(vj) and max(vi 3 ) > succ(vj). Consider the set of 
values Vi £ [v^^v^]. From Lemma fj\ Vj is the smallest 
available support for each of these Vi. Further, since Cij 
is a DS, max(vi) < succ(vj). Hence after pruning vj, no 
Vi has any support in Dj and is thus added to A. Also ob- 
serve that since max(vi 3 ) > succ(vj), any v[ £ [vi 3 , Vi 4 ) 
has succ{vj) as a valid support after pruning Vj. Hence, 
no v[ is added to A. □ 



5.2.5 Complexity of ACiDS Algorithm 

We first analyze the time complexity of ACiDS. The sub- 
routine Initialize() runs in 0(cd) time. For each constraint 
Cij £ C, this procedure simultaneously scans the domains 
Di and Dj without revisiting any value. Next, we note 
that ArcCons(i,j, A) runs in 0(d) time. From the de- 
scription of LocalArcConsQ in Section]^] we argue that 
it requires 0(A) time. 

It was stated in [14] that if the time complexity of Ar- 
cCons() is O(d) and that of LocalArcCons() is 0(A), 
then the algorithm AC5 runs in 0(cd) time. Since the 
ArcConsQ and LocalArcConsQ in ACiDS take O(d) 
and 0(A) time, we conclude that ACiDS also runs in 
0(cd) time. 

Analyzing the space complexity, we note that since 
each DS constraint can be stored as a list of intervals, the 
total input space requirement for a DSCSP is 0(cd). Dur- 
ing the execution, since any Vi £ Di is enqueued at most 
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Algorithm 3 Subroutines for Arc Consistency for USCSP Algorithm 4 DSCSP Solver 



procedure Initialize 

Vvj e Dj, VCy e C, MIN (y jt i)<- <j> 
for each Cy S C do 

scani , scanj point to last and first values of Dj , Dj re- 
spectively 

while scani ^(/>A scarij ^ do 
while min(scarii) > scarij do 

scarij <r- min(scam) 
MIN (scarij, «)<— Mi" TV (scarij, i) U {scan*} 
scam <— pred (scarii) 

procedure LocalArcCons(j, j, u; , A) 
if MIN(vj,i) = <£then 
return 

A^ </> 

MIN (vj ,i) = [t^ x , u, 2 ] for some Uj. , t>i 2 £ D, 

f' <— succ(vj) 

scarii 4— Vi„ 

while ( v'j = </> V moi (CVj , scan;) < v'. 1 A (scarii > "»i ) 



do 



A-s— A U {scarii} 
scarii <— pred(scarii) 
if scarii 7^- <P A scarii > ^h then 

M JiV (v< , i J <-MIN (v'j ,i) U[v tl , scarii] 



once for each arc(j, i) £ arc (V), the queue length of Q 
is at most 0(cd). By a similar argument, the MIN data 
structure takes at most O(cd) space. Since at any point, 
A holds values of any one domain, A takes no more than 
0{d) space. Hence, the ACiDS algorithm uses at most 
O(cd) space. 

The 0(cd) space and time complexity is optimal for 
any CSP with c constraints, since each constraint needs 
fl(d) space, and at least one scan of any domain requires 
fi(d) time. The ACiDS algorithm improves the complex- 
ity bound of 0{cd 2 ) in [36) for the subclass DS of CRC 
constraints. 

5.3 Arc Consistency for USCSP 

In this section, we show that with minor modifications 
in the way we scan the domains in the subroutines, Ini- 
tializef) and LocalArcConsf), the ACiDS algorithm can 
compute arc consistency for a USCSP. The modified sub- 
routines are presented in Algorithm [3] In addition to the 
assumptions stated in the ACiDS algorithm, we assume 
scarii = <fi if and only if scarii reaches the end of the do- 
main on either side, and pred(vi) = <fi if and only if Vi is 
the smallest value of the domain. It is easy to follow that 
arc consistency for a USCSP can also be achieved with the 
same time and space complexity of 0(cd) as in the case 
of DSCSP. However, since arc consistency is not known 
to be sufficient to determine a solution to a USCSP, tech- 
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procedure mainCT 3 = (X,V,C)) 
InitializeQ 

while Vi scarii 7^ <j> do 
j ^Dequeue (Q) 
for each {k : Cj k S C} do 

while scarij 7^ <j> A scarij ^ Dj (Cjh) do 
flagj = now 

Increment scarij to the next value in Dj 

while scan k 7^ <j> A scarij ^ D k (C k j) do 
flag k = now 

Increment scarij to the next value in D k 
if scarij < rain (C k j, scan^ then 

scanj <— min (C k j, scan k ) 

continue 
else if scan j. < min (Cjk, scanj 

scank <— min (Cjk, scanj) 

if flag k = 1 then 
flag k «- 



\> flagj = now 



t> Restart for loop 
then 



1 



flagj 



Q<-Qu{fc} 
1 



sum <— sum + 1 
if sum = n then 

print Solution = {Vi, scani} 
break 
procedure Initialize 
sum <— 
Q<- 
for each Xi £ X do 

scani points to the first value of D; 
flagi <r- 
Q ^ QU{i} 



niques such as path consistency or variable elimination 
are required. Since arc consistency is a basic subroutine 
in these techniques, replacing that by our modified ACiDS 
algorithm produces a faster solution. It improves the run- 
ning time of O {a 2 nd + cd 2 } (where a < n is a problem 
specific parameter) of the best known algorithm |[36l to 
O (a 2 rid + cd), which is linear in d. 



6 A Faster Solver for DSCSP 

This section presents a more efficient algorithm DSCSP 
Solver (Algorithm 0} to solve a DSCSP than ACiDS. 
We observe that a solution to a DSCSP can be found 
even without achieving arc consistency. Exploiting the DS 
property of the DS constraints and the total ordering of 
the variable domains, the algorithm incrementally scans 
the variable domains in a manner such that no value is re- 
visited. At each point, we check whether the values being 
scanned form a solution. Only when it is confirmed that 
the current value cannot form a solution is the next value 
scanned. 
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The algorithm maintains a queue Q of variable indices. 
A variable index i £ Q if and only if it is needed to check 
whether (scani, scarij) £ Cij for all the constraints Cy, 
imposed on X t . The algorithm also uses pointers scarii 
and boolean variables flagi for each variable Xi £ X. 
The variable sum maintains the sum of all the flag vari- 
ables. As stated earlier, we assume that for each pair of 
variables (Xi,Xj), there is at most one constraint Cy, 
and we need not consider its transpose Cji separately. (In 
the ACiDS algorithm, we considered both C^, Cji £ C.) 

DSCSP Solver first calls Initialize() for initializing the 
scarii pointer to point to the first value of Di and flagi to 
for each Xi £ X. The queue Q is initialized to contain 
all the variable indices. 

Then, the algorithm proceeds in iterations. In each it- 
eration, we dequeue a variable index j and consider all 
the constraints Cjk £ C. We keep incrementing scarij 
(and scarik) until it reaches a value which has a support 
in Dk (Dj respectively). If scarij < min(scarik), we set 
scarij = min(scarik) and reconsider all the constraints 
imposed on Xj. Also, if scarik < min(scanj), we set 
scarik = min(scarij). Further, if flagk = 1, we reset it 
to and add k to Q. Since sum maintains E"=i /^ a <?i> 
sum is decremented by 1 as a consequence. At the end 
of the for loop (lines 5|T8 1, if scarij , scarik ^ <t>, then 



0. Later, flagj is set to 1 in line 19 And immediately, 



sum = n, it implies Vj flagj = 1, which indicates 
that (scarij, scarik) £ Cjk for all constraints Cjk £ C. 
This in return signifies that at this point, the assignment 
Xi = scarii is a solution to the DSCSP. 



6.1 Correctness of DSCSP Solver 

To prove the correctness of the algorithm, we first show 
that the answer reported by the algorithm is indeed a so- 
lution. The following simple claims lead to this fact. 



Lemma 9. At the end of the for loop in lines 5\18 
flagj — 1 Vj if and only if j ^ Q. 

Proof. Consider a variable index j in Q. After j is de- 
queued from Q, the for loop in lines 5fT8 is processed. 
At the end of the for loop, flagj is set to 1, On the other 
hand, before any variable index k is added to Q in line 18 
flagk is decremented to in line 16 □ 



Lemma 10. At any stage of the algorithm, sum — 
E"=i f la 9i- 

Proof. From the previous lemma, it is clear that when any 
variable index j is dequeued from Q in line El flagj = 



sum is also incremented. On the other hand, when for any 
variable index k, whenever flagk is decremented from 1 
to in line 16 sum is also decremented by 1. □ 



Lemma 11. When flagj moves from to 1 (in line 191 



then (scarij, scarik) £ Cjk for each constraint Cjk on 



Xj. 



Proof. At a high level, the proof follows from the obser- 



vation that before flagj gets incremented in line 19 every 



constraint on Xj is verified to be true for the current value 
of the scan pointers. 

Consider the for loop in lines [5p8] Consider any con- 
straint Cjk £ C. Note that during the execution of the 
algorithm, any scan pointer is never decremented. Fol- 
lowing the down staircase property of Cjk, it is clear that 



if (scanj, scank) (£ Cjk such that , 



£ Dj(Cjk) and 



scank € Dk(Ckj), either of the following two cases are 
possible but not both: (i) scanj < min (Cf-j, scank), or 
(ii) scank < min (Cjk, scanj). 

From the algorithm, in case (i), scanj is incremented 
to min (Ckj, scank). This ensures that (scanj, scank) £ 
Cjk- However, since scanj has been modified, we need 
to reconsider all the constraints Cjk on Xj. Arguing simi- 



for scanj, every constraint Cjk imposed on Xj is sat- j mes 
isfied by scank- As a result, flagj moves from to 
1 and, therefore, sum is incremented by 1. Whenever 



larly for case (ii), it is clear that at the end of the for loop in 

5p8 when flagj is set to 1, (scanj, scank) £ Cjk 

for each constraint Cjk on Xj . □ 



Corollary 1. At any stage during the execution of the al- 
gorithm, (i) for any constraint C jk, if both flagj = 1 and 
flagk = 1> then (scanj, scank) £ Cjk, and(ii) if sum = 
n, then there is a solution, namely, Vi Xi — scani. 

Now we show that the algorithm would never miss a 
solution, if it exists. The key point to note here is that 
scani is not incremented until it is confirmed that the cur- 
rent value (being pointed by scarii) cannot be a member 
of any solution to the DSCSP. The following lemma for- 
mally proves this claim. 

Lemma 12. If Vi £ Di is a member of a solution, 
then at every stage during the execution of the algorithm, 
scani < Vi- 

Proof. Suppose Vi is the first value which got crossed by 
scani during the execution of the algorithm even when 
it is a member of a solution S. Suppose the solution is 
S = {Vfc, X k = v k }. 

Consider the time when scani is about to be incre- 
mented from Vi. Observe that at this time, scani = v. L 
and scanj < Vj for all j ^ i. The value v.i can get 
crossed if and only if for some constraint Cij, either 
Vi 4- Di(Cij) or Vi < min(scanj). If Vi $_ Di(Cij), 
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then Vi cannot be a member of a solution, and hence, this 
is a contradiction. Considering the latter case, suppose 
scarij = v'j < Vj. Since (vj,Vi) G Cji, min(vj) < Vi. 
However, min{v'A > v t > min(vj) contradicts the DS 
property of Cji. □ 

Hence, we have shown that the solution reported by the 
DSCSP Solver is valid and that it would never miss to 
detect a solution, if one exists. This concludes the cor- 
rectness analysis. 

6.2 Complexity of DSCSP Solver 

From the above results, it follows that if Vi,u, G Di is 
the smallest value in Di which is a member of a solution, 
then the algorithm would report the solution Vi, X^ = V{. 
The complexity analysis leverages this fact. 

Theorem 3. // Sj is the degree of node Vj (corre- 
sponding to variable Xj) in the constraint network, and 
v s is the smallest value in Dj that is a member of a 
solution, then the algorithm DSCSP Solver takes time 

O ( Sy=i &i s i ) ~ O ( cs ) w ^ere s — maxj {sj}. 

Proof. For a given variable Xj that is dequeued from Q, 
the maximum number of iterations of the for loop in lines 
|5|T8| is Sj as the number of increments of scarij is at 



most Sj — 1. Now consider any given iteration of this for 
loop for the variable Xj. Suppose scarij — Vj. Before 
scarij gets incremented during this iteration, Vj is com- 
pared with at most scarik for all k such that Cjk G C. 
This implies that in this iteration, Vj is compared with at 
most Sj values. Hence, the result follows. 



□ 



In the absence of any solution, the algorithm takes 
0(cd) time because at least one of the scan pointers must 
scan the entire domain. Analyzing the space complexity, 
besides the input which takes O(cd) space, this algorithm 
requires only 0(n) additional space for storing the scan 
pointers. 

Comparing the DSCSP Solver with the ACiDS algo- 
rithm, the key advantages are the following: (i) DSCSP 
Solver is much simpler to implement, (ii) it need not scan 
the entire domain if there is a solution, and hence, offers 
an improved time complexity of 0($^?=i SjS j ) as op- 
posed to O(cd) time in presence of a solution, and (iii) its 
additional space requirement is 0(n) as opposed to 0(cd) 
in the ACiDS algorithm. 

7 Conclusions 

This paper investigated the class of generalized staircase 
(GS) constraints, which generalizes the class of staircase 



constraints, while being a subclass of CRC constraints. 
We explored the properties of the two subclasses of GS, 
namely, down staircase (DS) and up staircase (US) con- 
straints, besides studying their relationships with the ex- 
isting constraint classes. 

We proposed two algorithms for solving DSCSP, 
namely, ACiDS, based on arc consistency, and DSCSP 
Solver, based on incremental scan of domains. We 
showed that DSCSP Solver is more efficient as compared 
to ACiDS. For solving USCSP, we require either path con- 
sistency or variable elimination technique. Our optimal 
ACiDS algorithm acting as a subroutine in these tech- 
niques, offer a faster solution to USCSP. 

In the future, more efficient arc consistency scheme for 
CRC constraints may be investigated. In addition, inter- 
esting generalizations of GS constraints and applications 
to real life problems will be explored. 
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